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Abstract 

We consider an ideal closed stock market, in which 100 traders have economic activities. The 
assets of the traders change through buying and selling stocks. We simulate the assets under 
conservation of both total currency and total number of stocks. If the traders are identical, then 
the assets are distributed as a stationary Gaussian. When variety among the traders makes winners 
and losers, the asset distribution displays power law scaling such as the Pareto law We discuss 
this power law scaling from the point of view of superstatistics P|. It is given as a superposition 
of scaled distributions for each hierarchical level. The various traders have the same growth rate 
distribution to keep the scaling. 
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More than 100 years ago, Pareto found the power law scahng in individual income 
distributions 1]. How does the universal statistical law appear from complicated economic 
activities? Many approaches have been proposed to answer this question Re- 
cently mean field dynamics with multiplicative noises explained the power law behavior 
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But individual characters were completely neglected in the mean field ap- 
proach. We reconsider the power law scaling from the point of view of superstatisticsO]. It 
is a result of the superposition of scaled typical statistical distributions with respect to the 
characteristics. 

We will introduce a model of the stock market to understand the behavior of asset 
distributions in this letter. The model is an ideal closed stock market. There is no flow 
in or out of traders, cash or stocks. The total number of traders, total cash and total 
number of stocks are kept constant. The rules of trade are given by the extended threshold 
modelQ. 
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In the comparison of finance and physics, we can find several common features 
In a stock market, a trader corresponds to a molecule in a gas system. Assets play the role 
of energy. The conservation of total assets means energy conservation. When traders are 
identical, we can expect the assets to become distributed as the Maxwell-Boltzmann. 

A different behavior is expected when the characters of traders are widely divergent. 
When there are superior and inferior traders in a stock market, some hierarchical levels such 
as poor, middle, and rich appear. The asset distribution displays power-law scaling, which 
is different from the Maxwell-Boltzmann type distribution seen for identical traders. As 
will be shown later, the power law scaling is realized as a result of superposition of scaled 
distributions for each hierarchical level. 

In this letter we define superstatistics as a superposition of scaled statistical distributions 
with respect to a characteristic parameter. Our definition may be more general than the 



original one by Beck-Co^ 
developed in some fields 



len 2L The i 

3,y,Q. I 



dea of a superposition of scaled distributions has been 
Recently Beck and Cohen introduced superstatistics to 



explain the Tsallis statistics given as a superposition of the Boltzmann distributions with 
respect to certain temperature distribution j^,]- They also consider more general statistics for 
another temperature distributions. 

From the aspect of superstatistics, the appearance of hierarchical structure is a key issue. 
A trader stays at a fixed hierarchical level. The possibility to be rich or poor beyond that 
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hierarchical level is quite small. This heterogeneous (non-ergodic) behavior is an important 
difference from the mean field approach. We will discuss the relation between the two aspects 
in the end. 

Now we consider an ideal closed stock market, in which traders buy and sell stocks. For 
simplicity, we choose only one brand of stock. The ith trader has currency Ei{t) and stocks 
Vi{t) at time t. The stock market is completely isolated. There is no flow in or out of traders, 
currency or stocks. Therefore, the number of traders A^, total currency Etot = '^f=iEi(t) 
and total number of stocks Vtot = J2iLi ^«(^) ^^P^ constant. 

From a statistical mechanical point of view, the assets of traders play the role of energy 
in the stock market. The ordinary assets of the ith trader is defined as Xj ~ Ei{t) + Vi(t)p{t) 
where p{t) is the price of stock at time t. In that definition, the total assets depend on time. 
In this letter we choose the following definition: 

Xi = E,{t) + V,{t)p{0). (1) 

Although this definition is unusual, we give priority to strict conservation of total assets 

Q 

The asset of the ith trader changes through buying and selling stocks. In a trade at time 

t, 

E,{t) = E,{t - 1) ± f,,{t)p{t) 
V,{t) = Vi{t-t)T|^^{t) 

where fii{t) is the volume of trade for the ith trader. The buyer/seller pays/receives currency 
fii{t)p{t) to buy/sell fii{t) stocks. Although the assets of the traders change, the total assets 
are kept constant. 

Each trader demands/offers a bid/ask price and a buy/sell order volume at every time- 
step. The bid and ask prices follow the rules introduced by Sato-Takavasu jl7 1 . The Sato- 
Takayasu model is a deterministic threshold model 3]. Although the model is very simple, 



its stock price has similar statistical properties as real stock prices 
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The bid/ask price of the ith trader is determined by the following rules: 

(3) 



Bi{t + 1) = Bi{t) + ai{t) + Ci\p{t) - p{tp,,,) 



Si{t + 1) = B,{t + 1) + Ai 
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where Aj, \ai\ and q are characteristic parameters of the ith trader. The bid/ask price 
depends on the up/down movement of the stock price Cj[p(t) — p(tprev)] where p(tprev) is the 
stock price of the previous trade. 

The sign of aj(t) changes as — aj(t — 1), when ai{t — 1) > for the buyer or when 
aj(t — 1) < for the seller. Otherwise the sign is the same. The bid/ask price changes 
in each step. After the ith trader becomes a buyer/seller, reverses its sign. As a result, 
the ith trader alternately becomes a buyer and seller. The period is approximately given as 
Ti ^ 2Ai/\ai{t)\. 

The Sato-Takayasu model misses the volume of trade, which we need to estimate the 

n 

asset distribution. We extended the threshold model to introduce a volume of trade 
We assume that a buy/sell order volume is proportional to the trader's assets. 

fi^iB,) = KE,it-l)/p{t~l) 

(4) 

/is(S'i) = KVi{t - 1) 

where fih{Bi) / fis{Si) is a buy/sell order volume with the bid/ask price of the ith trader, k is 
a constant. We choose < k, < 1. This means that no trader can demand/offer a buy/sell 
order volume beyond their own assets. 

We determine the stock price and volume of trade using the following cumulative demand 



and offer functions 
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Oiq) = E5,<g/^s(5i) 

The balance between the demand and offer determines the stock price and volume of 
trade. (1) If D{p*) = 0{p*) for a unique p*, then the stock price is p{t) = p*. If D{p*) = 
0{p*) is unique, the volume of trade is D{p*) = 0{p*). If D{p*) = 0{p*) is not unique, 
the volume of trade is the maximum of D{p*) = 0{p*). (2)If D{p*) = 0{p*) for all 
p* in [pminjPmax], then the stock price is p(t) = (pmin +Pmax)/2. The volume of trade is 
D{p*) = 0{p*). To satisfy the balance with respect to volume, we need to reduce some 
buy/sell order volumes. 

The change of the bid/ask price of the ith trader, |aj|, chracterizes the relative merit of 
the ith trader in the extended threshold model, since a trader with a smaller change can 



buy/sell stocks at a relatively lower/higher price. We will discuss this property in 22 1. 
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FIG. 1: Cumulative distribution of assets for common r. The dashed hne is the error function. 
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FIG. 2: Gumulative distribution of assets for exponentially distributed r. The dashed line is 

We simulated the time evolution of the stock market and obtained the asset distribution 
for long times. We chose the following parameters and initial values: = 100, Aj = A = 
1000, po = 3000, Vtot = 10000, E^^^ = po14ot, = 0.5, q = 0. The initial bid price of the ith 
trader Sj(0) is given as a random number within (po — A,po)- The initial ask price of the 
ith trader is given as 5'i(0) = -Bi(O) + Aj. Each trader has equal initial assets. 

We considered the following two typical cases with respect to the period of trade, Tj = 
2A/|a,| 

(i) Common Ti 

In this case, the traders are completely identical. The period of trade is almost the same. 
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FIG. 3: Average, standard deviation and higher order moments have the same exponent with 
respect to asset size. The points for each moment are distinguished according to their asset size 
into twenty groups to clarify the size dependence. 

There are no strong/weak traders. The mechanism of stabihty works in this financial market. 
A trader who orders a greater/less volume has a tendency to lose/gain the asset. Each trader 
gains and loses assets under the mechanism of stability. We numerically confirmed the assets 
to be distributed as a stationary Gaussian illustrated in Fig^ 
(ii) Exponentially distributed Xj 

We chose Xj as an exponentially distributed random number. In this case, winners and 
losers appear. Although the distribution function cannot be completely steady, it displayed 
the power law scaling illustrated in Figj^js^. 

The assets were distributed as a stationary Gaussian for identical traders. For diverse 
traders, the asset distribution displayed the power law scaling. These results of our numerical 
simulations suggest that some hierarchical levels such as poor, middle, and rich appear. And 
the power law scaling is given as a superposition of scaled distributions for each hierarchical 
level. 

We will explain how the power law appears from superstatistics for simple situations. 
The relative merit of traders is characterized by the period of trade = 2A/|aj|. Traders 
who have longer period become winners. For almost steady state, the time average of the 
assets of the ith trader is scaled as rrii ~ mexp[Arj]. Furthermore the standard deviation of 
the time average of the assets is similarly scaled as cTj ~ o"exp[Arj]. This scaling may come 
from the fact that offer volume is always proportional to a trader's asset. 
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r = Log[X(t)/X(t-40)] 

FIG. 4: Probability density of growth rate r = log(X {t)/X{t — 40)). The points are distinguished 
according to their asset size into five groups. 

If the time average of the assets of each trader is distributed as a Gaussian, the time 
averaged asset density P{x) can be derived as 



where /(t^) ~ e is the probability density of the ith trader. It is not necessary a Gaussian. 
Any scaled distributions are sufficient. 

A= . Ref.a. both .a,ed di*ibu«on and detailed balance a.e efficient 

conditions for the Pareto law. To confirm our approach we calculated the mean, the standard 
deviation and higher order moments. These quantities showed scaling behavior as is shown 
in FigOl The Gibrat law (the size independence of the growth rate) is also illustrated in 
FiglU Notice that both the winners and losers have the same growth rate distribution to 
keep the Pareto law, even they have the different abilities. These results are consistent with 
the theorem 0, y,|3- 

In this letter, we have considered a model of the stock market in which 100 traders have 
economic activities. The assets of the traders change through buying and selling stocks. We 
have simulated the assets under the conservation of both total currency and total number of 
stocks. The assets become distributed as a stationary Gaussian for identical traders. When 
differences among traders make winners and losers, the asset distribution displays power law 
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FIG. 5: The large triangle with □ and x represents the region in which Eqs.(7) are satisfied. If 
\ai\ < \aj\, the area of Ap{t) > (x) is less than Ap{t) < (□). 

scaling such as the Pareto law. 

We have interpreted this power law scaling from the point of view of superstatistics. 
When the diversity of traders makes superiors and inferiors in a stock market, some hier- 
archical levels such as poor, middle, and rich appear. The asset distribution is given as a 
superposition of scaled distributions for each hierarchical level. 

From the viewpoint of superstatistics, the appearance of hierarchical structure is a key 
issue. A trader stays at a fixed level. The possibility to be rich or poor beyond the cor- 
responding level is quite small. This heterogeneous (non-ergodic) behavior is an important 
difference from the approach of mean field dynamics. 

For an open market, the abilities of traders change in time. Strong traders become weak 
and vice versa. Only traders who can fit the present market become winners. Therefore the 
hierarchical structure becomes fuzzy for long time scales. We expect that the mean field 
dynamics are recovered for such long time scales. 

For income distributions of companies, both the Pareto and Gibrat laws are valid in some 
well developed countries such as USA. and Japan . As we have showed, the size 

independence of the growth rate does not deny the existence of hierarchical structures. The 
probability to be above or below the corresponding size is much smaller than what the mean 
field dynamics predicted. The hierarchical structures can be confirmed from empirical time 
series data. Growth rates can be related to the abilities of companies in quickly developing 
countries. 
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[22] Suppose the jth trader sells some stocks and the ith. trader buys them at time t. The ask price 
of the jth trader is lower than the bid price of the ith. trader: Bi(t) = Bi(t — l) + \ai\ > Sj{t) = 
Sj{t — 1) — \aj\. Furthermore the stock price is lower than the ask price of the j'th trader and 
lower than the bid price of ith trader at time i — 1: Bi(t — 1) <p{t — l) < — 1). From these 
relations, we obtain ABi{t - 1) < 0, ASj{t - 1) > 0, and ASj{t - 1) > ABi{t - 1) + |ai| + \aj\ 
where ABi{t -I) = Bi{t - 1) -p{t - 1) and ASj{t - 1) = Sj{t - 1) -p{t - 1). As is illustrated 
in Fig|SJ the region in which the above three conditions are satisfied is given as a triangle. We 
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If |aj| < \aj\, the area of Ap{t) > is less than Ap(t) < in FigO Therefore the ith trader 
can buy stocks relatively lower price. If |aj| > \aj\, the jth trader can sell stocks relatively 
higher price. 

[23] Although we discussed the power law scaling, the asset distribution is not completely steady. 
This comes from the instability that rich traders have a tendency to gain big benefits even 
under the conservation of total assets. The picture of equilibrium with small fluctuation is 
not valid. There always exist large fluctuations. To stabilize the distribution the role of tax 
becomes important. We take some cash from rich traders and distribute them to poor traders. 
We consider that a good tax keeps the scaling of distribution, and a bad one breaks it. 
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